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Abstract 

The Ising spin - S model on recursive p - polygonal structures in the external magnetic 
field is considered and the general form of the free energy and magnetization for arbitrary 
spin is derived. The exact relation between the free energies on infinite entire tree and on 
its infinite "interior" is obtained. 
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1. INTRODUCTION 



Many theoretical problems of statistical mechanics, solid state physics, gauge models, poly- 
mers, etc., can be solved exactly on recursive structures, like the Bethe |], H, |3|, HI an d the 
Husimi |5|, ^ |7J lattices. Moreover, it has been argued recently || that, in some case, the 
Bethe lattice calculations are more reliable than mean-field calculations. In connection with 
this, the knowledge of the exact form of the free energy and magnetization is very actual for 
locating phase transitions and for analytical investigations of physical phenomenon. The main 
problem is how to calculate the free energy on the infinite "interior" of the infinite entire lattice. 
Namely, how get rid of the surface (with nonzero density of sites) contributions which can allow 
to strange results and rather different from those on regular lattices 0. In our previous papers 
the relation between the free energies on Bethe and Cayley trees for spin-1 Ising model || and 
on Husimi and corresponded to it Cayley trees for multisite Ising model J7| are obtained. In 
that approach the free energy derived from the recursion relations and then the differentiation 
of the free energy functionals over the external magnetic field allowed to the exact relations. 
Recently, a method of computing the bulk free energy for any model on infinite "interior" tree 
is presented @. 

In this paper we have generalized the approach of the surface independent free energy 
calculation [|3], [7| for an Ising spin - S model in an external magnetic field on generalized 
recursive lattice consisted of p - polygons (p is the number of edges of the polygon). 

2. THE MODEL 



The recursive tree consisted of p - polygons is characterized by p, the number of edges (the 
number of sites) of the polygon (p = 2 - usual Bethe lattice, p = 3 - Husimi lattice) and by q, 
the number of p-polygons which go out from each site (Fig.l). 

These infinite dimensional recursive lattices cannot be embedded in any lattice belonging 
to an Euclidean space of finite dimension. However, they can be embedded in two-dimensional 
space of constant negative curvature (the hyperbolic or Lobachevsky plan H2) and study their 



metric properties 10 



Let us define the Ising spin-S" model in an external magnetic field by Hamiltonian 

-PH= £ H l {a ( t\af\---af ) )+Y, H Ml (1) 

p—poligon i 

where cr, take values (— S, —S + 1, S — 1,S), - spin around polygon, the first sum goes 
over all p - polygons and the second over all sites.. H\{a^ \ af \ a\ p ^) - includes all possible 
nearest - neighbor pair interactions, as well multisite interactions between spins belongs the 
same polygon and H 2 ((Ji) = Y^fLi hp. a i ~ includes all possible single - ion interactions. 
The partition function will have the form 

■■^n • z^n, (2) 
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The advantage of the recursive tree is that for models formulated on it exact recursion 
relations can be derived. When this tree is cut apart at the central site 0, it separates into q 
identical branches, each of which contains (q — l)(p— 1) branches. Then the partition function 
can be written as follows: 

Z = ^exp{H 2 (<To)}g%(<T ), (3) 

CO 

where <jq is a spin in central site, N is the number of generations, and gjv(cr ) is the partition 
function of a branch. Each branch, in turn, can be cut along any site of the lth-generation 
which is nearest to the central site. The expression for g^{ao) can therefore be written in the 
form 

9n(*o) = E expllfiK^...*^ (4) 
{*?>} I i=1 ) 

where are the spins of first polygon except the central site. Consequently, we will have 
2S + 1 recursion relations for the gN^o), where a take values (—S,—S + 1,....,S — 1,S). 
After dividing each of the recursion relations on a recursion relation for g^(S), we will have 25* 
recursion relations for xn(cto) 



x N {a 



£ {CT o> } exp {i^K a« • • • a^) + E^l H 2 {af)} x%l x {a?) ■ ■ ■ x^a 



N-1K U 1 



E {CT u) } exp faS, a?\ ■ ■ ■ a™) + Eg H 2 (a?)} x^a?) ■ ■ ■ 

where 



N-l\ a l ) 



(5) 

Xn(ct ) = g N (a )/g N (S), (6) 

and an equation for Qn(S) 

gN(S)=g i r-l ]ip - 1 \S)^N-i, (7) 



where 



*s-i= E e^\H 1 (S,*[ 1 \---*t 1) )+ P j:M*? ) )\^\(^ (8) 

Since the right-hand side of Eq. (g) is bounded in x^, it follows that x^ is finite for iV — > <x>. 

Through this xat(ct) one can express the density m M = (<Jq) of central site (the symbol (...) 
denotes the thermal average), where [i takes values from 1 to 25": 

m _ _ Ea go exp {g 2 (a )} x^(a ) 

''V — \ a 0/ — 7TT~i vi <T7 \~ \ J ) 

and other thermodynamic parameters. So we can say that the xn{o) in the thermodynamic 
limit (N — > oo) determine the states of the system. 

Using the Eq. (|(J), the free energy on the Cayley tree can be written as follows 

- f3FS ayley = q\ng N (S) + \n($ N ), (10) 
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where 

^ N = ^2expH 2 {(To)x%(<TQ). (11) 

CO 

By substituting in the right-hand side of Eq. (|TUD the expression of gjv(S') given by Eq. (|7J), 
we obtain the recursion relation for the free energy on the generalized Cayley tree 

- [3Fg ayley = -(q - l)(p - VPFfft*" + qH**-i) ~ (ff " 1)(P " 1) ln(*«r-i) + (12) 
By repeating this recursion procedure n times one can obtain 

- (3FS aylev = -[{q- l)(p - ^ - /3F nJV , (13) 

where 

n 

- /3F n7V = q £ [(? - 1)(P - i)]"" 1 In(^iv-ic) - [(? - l)(p - l)] n In(^-n) + ln($iv). (14) 
fc=i 

Therefore we obtain the general form of the free energy for Ising spin-S* model on generalized 
Cayley tree 

N-l 

- (3FS ayley = q £ [(? - 1)(P - I)]*"*- 1 H**) + In(^) (15) 

fc=0 

where ^ and &n are given by Eqs. @ and (|TT|) . 

As mentioned in the Introduction the above obtained free energy gives rise to unusual 
behavior due to surface sites (in particular, it depends on surface conditions) 0. To overcome 
this problem we here consider only local properties of sites deep within the graph (i.e. infinitely 
far from the boundary in the limit N — > oo). Indeed, we will consider only contribution to 
from the Bethe lattice. 

Let us consider the case, when the series of solution of recursion relations given by Eq. (|5|) 
converge to a stable point at N —>■ oo. In this case limxjv_ n (cr) = x(cr) for all finite n and the 
recursion equations become 

E {g o) } exp jffxK gg, • • ■ + Eg H 2 {a^)} ^-\a?) ■ ■ ■ x^ja^) 

x{ao) " z^^i^^ ' (16) 

Then the Eq. (|14}) will take the form 

-(3F n = lim {-(3F nN ) 

N^oo 

(a - l) n (v - l) n - 1 
where \& = liniTv^oo ^N-k and $ = liniAr^oo &N-k for all finite k and 

* = £ exp <r?\...<^^ ( 18 ) 



$ = 5>xp# 2 (a )xVo). (19) 

Let us now prove that Eq. fll7l) gives the free energy of the spin - S Ising model on the 
generalized Bethe lattice consist of n generations. 

First, let calculate the total number of the sites of the generalized Bethe lattice consist of 
n generations. By definition, the Bethe lattice is the union of equivalent sites and therefore it 
is easy to get the relation between the numbers of bonds and sites. As from each site go out q 
bonds and each bond contains two sites then N bond = ^N site . On the other hand, the number 
of bonds of the Cayley tree or the Bethe lattice consist of n generations is 

n-i (Q -l) n -l 



N bond = q + q(q-l) + q(q-iy + ■■■ + q(q - l)™" 1 = q 



q-2 

So, the total number of sites of the Bethe lattice consist of n generations is 

AT _ 2 N -O P?' 1 )"- 1 

1 v site 1 > bond " n 

q q-2 

Above we calculate the the number of sites in the case when p = 2. It is easy to generalize this 
result for arbitrary p. Note, only, that in case when p > 2 from each site go out 2q bonds and 
the relation between the numbers of bonds and sites will looks like 

2q 

Nbond = —N s ite = qN site , 

and the number of bonds for arbitrary p will have the following form 
N bond = qp + qp(q - l)(p - 1) + • • • + qp(q - l) n ~\p - l) n_1 = qp 



(g-l)(p-l)-l ' 

Thus, the total number of sites of the generalized Bethe lattice [] consist of n generations is 

]\[Bethe _ (g - l)"(p ~ 1)" - 1 { ) 

IVstte ~ P (g_l)(p-l)-l ' (M) 

Now we can write the explicit expression for the free energy per site for the spin - S Ising 
model on the generalized Bethe lattice 

- = -^g, = q - M») - {q ~ - 1} - 1 ■■(•)■ (2i) 

ly site P P 



1 One can reach to the same result in a different way as well. Indeed, one of the ways to fulfill the equivalent 
of the sites of the Bethe lattice is to put a periodical boundary condition on Cayley tree. Then the number of 
sites of the last shell of Cayley tree will be equal to q(q — — l) n_1 , whereas the number of sites of the 

corresponding Bethe lattice is (q— l) n ^ 1 (p— 1)™ . One can see that the number of sites on the boundary shall 
of the generalized Bethe lattice is q times less then the corresponding number on the generalized Cayley lattice. 
Consequently 

N*tl he = 1 + q(p - 1) + q(q - l)(p - l) 2 + • ■ ■ + q(q - l)»- 2 (p - l)"" 1 

(«-l) B (p-l) n -l 



-{q-ir-\p-ir= P - 



(,-i)(p-i)-i 
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where \& and <£> are given by Eqs. ([Tj|) and (|T9|) respectively. 

We wish to mention here that in our earlier paper we had noted that this functional 
gives the exact form of the free energy for spin- 1/2 Ising model on the Bethe lattice (p=2) (Tl| . 



Then we checked the correctness of it for a spin-1 Ising model on the Bethe lattice (p=2) by 
differentiation over the external magnetic field. Recently, the correctness of that functional is 
confirmed for multisite antiferomagnetic spin-1 /2 Ising model on the Husimi tree (p=3) as well 

0- 

Below an effort is made to prove that this statement is correct for arbitrary spin - S Ising 
model on the generalized Bethe lattice. 

With this aim in view, let us differentiate this free energy functional with respect to the 
external field h^. Without loss of generality, we will consider, for simplicity, the case p = 2. 
For arbitrary p, the proof will be accomplished in the same way. 

After a little algebra, we obtained 

S , R f \ E ffo <rftexp{# 2 (cro)}xVo) 

-{P J Bethe j 



Here 



, g Ea exp{H 2 (a )}x^\a )T(a ) 

2E^ ex P {H 1 (S,a) + H 2 (a) + H 2 (a )}x^(a)x^(a y { > 



T(ao) = Eexp{^ 2 (^)K- 1 (a) 



x'(a) 



x l[a k + (q - 1)-^] [x(a ) exp {H^S, a)} - exp {H^q, a)}] + x'(a ) exp {H^S, a)} , (23) 



where 



It is easy to check that 



T((x ) = ^- S^£ex V {H 2 (a)}x^\a) [x(a )exp{H 1 (S,a)} - exp {i^K a)}] j . (24) 

After substituting in the right-hand side of Eq. (p4"|) the expression of x(cxo) (Eq. (|16"|)), which 
in the case p = 2 is 

( v _ Y, a ^{Hi{a^a) + H 2 {a)}x^ l {a) 
X[CT0) J2 fT exp{H 1 (S,a) + H 2 (a)}x^(a) ' 



we will obtain 
and consequently 



T(<7 ) = 



6 -fc) = m,= %^g^^g^ . (26) 



(tq ex P {-^(oo)} ^ 9 (c"o) 
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Thus, the expression given by Eq. ([H]) is the exact free energy functional per site for an 
Ising spin-S 1 model on the generalized Bethe lattice. Consequently, the equation (Eq. flT3p) Q is 
the exact relation between the free energies on generalized Cayley tree and generalized Bethe 
lattice for a spin - S Ising model. It is obvious that this relation is model independent and 
depends only on the structure of recursive tree. 

3. CONCLUSION 

In this paper we have considered an Ising spin - 5* model on recursive p - polygonal structures 
in the external magnetic field and derive the general form of the free energy and magnetization 
for arbitrary spin. We have obtained the exact relation between the free energies on infinite 
entire tree and on its infinite "interior" and show that it is model independent and depends 
only on structure of recursive tree. The advantage of approach presented in this paper is that it 
gives possibility to obtain the exact expression for the surface independent free energy for width 
class of spin systems on generalized recursive lattices. This approach should be applicable for 
gauge models on generalized multi-plaquette recursive structures as well.. 
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2 More recently, an elegant geometrical interpretation of this relation is presented in Ref. ||. The key idea 
of that method is the construction of a union trees with coordination number q from initial tree with the same 
coordination number. 
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5. Figure Caption 



Fig.l The recursive tree with p = 4 and q — 2. The numbers 0, 1, 2 and 3 denote the central 
site, the sites of 1th, 2nd and 3rd generations, respectively. 
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